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Abstract 

An effective hadronic lagrangian consistent with the symmetries of quan- 
tum chromodynamics and intended for apphcations to finite-density systems 
is constructed. The degrees of freedom are (valence) nucleons, pions, and the 
low-lying non-Goldstone bosons, which account for the intermediate-range 
nucleon-nucleon interactions and conveniently describe the nonvanishing ex- 
pectation values of nucleon bilinears. Chiral symmetry is realized nonlinearly, 
with a light scalar meson included as a chiral singlet to describe the mid-range 
nucleon-nucleon attraction. The low-energy electromagnetic structure of the 
nucleon is described within the theory using vector-meson dominance, so that 
external form factors are not needed. The effective lagrangian is expanded 
in powers of the fields and their derivatives, with the terms organized using 
Georgi's "naive dimensional analysis". Results are presented for finite nuclei 
and nuclear matter at one-baryon-loop order, using the single- nucleon struc- 
ture determined within the model. Parameters obtained from fits to nuclear 
properties show that naive dimensional analysis is a useful principle and that 
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a truncation of the effective lagrangian at the first few powers of the fields 
and their derivatives is justified. 



I. INTRODUCTION 

Quantum chromodynamics (QCD) is generally accepted as the underlying theory of the 
strong interaction. At low energies relevant to nuclear physics, however, nucleons and mesons 
are convenient and efficient degrees of freedom. In particular, relativistic field theories of 
hadrons, called quantum hadrodynamics (QHD), have been quite successful in describing 
the bulk and single-particle properties of nuclei and nuclear matter in the mean-field and 
Dirac-Brueckner-Hartree-Fock approximations [0,0 . 

QHD studies based on renormalizable models, however, have encountered difficulties due 
to large effects from loop integrals that incorporate the dynamics of the quantum vacuum 
On the other hand, the "modern" approach to renormalization P-[T(][], which makes 



sense of effective, "cutoff" theorie^ with low-energy, composite degrees of freedom, provides 
an alternative. When composite degrees of freedom are used, the structure of the particles is 
described with increasing detail by including more and more nonrenormalizable interactions 
in a derivative expansion §]. We seek to merge the successful phenomenology of the QHD 
approach to nuclei with the modern viewpoint of effective field theories. 

The present work is an extension of Refs. [|ll|-13|, which explored some aspects of incor- 



porating QCD symmetries into effective hadronic models. Here we build a more complete 
effective lagrangian and apply it to nuclei at the one-baryon-loop level. Our two main goals 
are to include effects of compositeness, such as single-nucleon electromagnetic structure, in a 
unified, self-consistent framework and to test whether "naive dimensional analysis" (NDA) 
and naturalness (see below) are useful concepts for effective field theories of nuclei. Eventu- 
ally, we will need to demonstrate a systematic expansion scheme that includes corrections 
beyond one-loop order. Nevertheless, Ref. [0 describes how a general energy functional 
with classical meson fields can implicitly incorporate higher-order contributions. We there- 
fore expect that higher-order many-body effects can be absorbed into the coefficients of our 
one-loop energy functional when we determine parameters by fitting to nuclear observables. 
We return to this point later in this section. 

Important physical constraints on the effective lagrangian are established by maintaining 
the symmetries of QCD. These include Lorentz invariance, parity invariance, electromagnetic 
gauge invariance, isospin, and chiral symmetry. These symmetries constrain the lagrangian 
most directly by restricting the form of possible interaction terms. In some cases, the 
magnitude of (or relations between) parameters may also be restricted. 

Chiral symmetry has been known to play an important role in hadronic physics at low en- 
ergies since the early days of current algebra |]14|,|15| ; it is exploited systematically in terms of 
effective fields in the framework known as chiral perturbation theory (ChPT) [|TB-|T8|. There 
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is also a long history of attempts to unite relativistic mean-field phenomenology based on 
hadrons with manifest chiral symmetry. In particular, it has been tempting to build upon 
the linear sigma model p!9|-pT|, with a light scalar meson playing a dual role as the chi- 
ral partner of the pion and the mediator of the intermediate-range nucleon-nucleon (NN) 
attraction. The generic failure of this type of model has been shown in Refs. [|ri|,|T^.0 How- 
ever, a light scalar is completely consistent with a nonlinear realization of chiral symmetry 
|]T^,|TB|. Moreover, in a nonlinear realization, one need not assume that the hadrons are 
grouped into particular chiral multiplets p4 |. 

The phenomenological success of ChPT establishes the dominant contributions of the 
vector resonances to the low-energy constants in the pion sector p5|-p^. In particular, 
a tree-level effective lagrangian with vector mesons is found to be essentially equivalent 
to ChPT at one-meson-loop order. In the baryon sector, a satisfactory description of the 
nucleon form factors from vector-meson dominance (VMD) [p8H 30[| alone has not been found 
Pl| . However, we can describe the long-range parts of the single- nucleon form factors (that 
is, magnetic moments and mean-square radii) using an effective lagrangian that includes 
vector mesons and direct couplings of photons to nucleons. Our approach is analogous to 
that implemented long ago in Ref. |P2|| . 

In virtually all previous QHD calculations of finite nuclei, the single-nucleon electromag- 
netic structure was incorporated "by hand." After the nuclear structure calculation was 
completed and "point-nucleon" densities were obtained, they were convoluted with empir- 
ical single-nucleon form factors to compute the charge density. The results were typically 
insensitive to the nature of the form factor, as long as the nucleon charge radius was ac- 
curate. In the present work, we do not need external form factors. The large oscillations 
in the point-nucleon charge density of a nucleus are significantly reduced by virtual mesons 
inherent in the model, and the low-momentum parts of the charge form factors of various 
nuclei are well reproduced. Thus we have a unified framework. 

To achieve a systematic framework for calculating nuclei, however, we must develop an 
expansion scheme. ChPT is an expansion in powers of small external momenta and the 
pion mass. Weinberg |^ applies this small-momentum expansion further to a calculation 
of the iV-nucleon "effective potential" , defined as the sum of all connected diagrams for the 
S matrix (in time-ordered perturbation theory) without A^-nucleon intermediate states. On 
the other hand, an expansion of the NN interaction kernel in terms of intermediate states of 
successively shorter range provides a successful and economical way to accurately describe 
NN scattering data |^|. Here we combine the ideas of the small- momentum expansion and 



the intermediate-state expansion by including the low-lying non-Goldstone bosons into the 
low-energy effective lagrangian. These bosons account for the intermediate-range interac- 
tions and conveniently describe the nonvanishing expectation values of nuclear bilinears such 
as NN and N'jf.N. 



^The failure of models based on the linear sigma model lies more with the assumed "Mexican-hat" 
potential than with the linear realization of chiral symmetry. Successful mean-field models with 
a linear realization and a light scalar meson have been constructed using a logarithmic potential 
(see Refs. ||,||). 
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There are well-defined resonances in the vector channels, and thus it is natural to identify 
them with fields in the effective lagrangian; in contrast, the dynamics in the scalar channel is 
more complicated. Numerous calculations p5|-p9| show that correlated two-pion exchange 
leads to a distribution of strength in the isoscalar-scalar channel at center-of-mass pion 
energies of roughly 500 to 600 MeV, even in the absence of a low-energy resonance. We 
therefore introduce an isoscalar-scalar meson as an effective degree of freedom to simulate 
this strong interaction, just like the non-Goldstone bosons in the vector channels. Such a 
scalar, with a mass of 500 to 600 MeV, is essential in one-boson-exchange models of the NN 
potential and in relativistic mean-field models of nuclei We emphasize that at the 
level of approximation used in this paper, the scalar field can be introduced without loss of 
generality. (For example, it can be viewed as a way to remove terms with more than two 
nucleon fields from a general pion-nucleon lagrangian.) A deeper understanding of the NN 
interaction in the scalar channel is an important pursuit, but is not needed for our purposes 
here. 

We shall expand our effective hadronic lagrangian in powers of the fields and their deriva- 
tives, and organize the terms using "naive dimensional analysis" (NDA) The NDA 
framework is used to identify the dimensional factors of any term. After these dimensional 
factors and some appropriate counting factors are extracted, the remaining dimensionless 
coefficients are all assumed to be of order unity. This is the so-called naturalness assumption. 
If naturalness is valid, the effective lagrangian can be truncated at a given order with a rea- 
sonable bound on the truncation error for physical observables. Thus we have a controlled 
expansion, at least at the tree level. 

Evidently, for an expansion in powers of the fields and their derivatives to be useful, 
the naturalness assumption must hold. Implicit in this assumption is that all the short- 
distance physics is incorporated into the coefficients of the effective lagrangian; otherwise, 
the truncation error will not be reliably given by NDA power counting when loops are 
included. Thus, just as in ChPT or in heavy-baryon theories, the nucleon vacuum loops and 
non-Goldstone-boson loops will be "integrated out", with their effects buried implicitly in 
the coefficients.^ Nevertheless, the nucleon and non-Goldstone-boson fields are still needed 
in the lagrangian to account for the valence nucleons, to treat the mean fields conveniently, 
and to apply the usual many-body and field-theoretic techniques. Any contributions from 
vacuum nucleon or non-Goldstone-boson loops that arise in many-body calculations can 
formally be cancelled by the appropriate counterterms. These cancellations can always 
be achieved, since all possible interaction terms consistent with the symmetries are already 
included in the effective lagrangian. On the other hand, long-distance, finite-density effects 
should be calculated explicitly in a systematic application of the effective lagrangian. 

In principle, the naturalness assumption can be verified by deriving the effective la- 
grangian from QCD. Since this is not yet feasible, one must fit the unknown parameters 



^We note that the one-nucleon-loop vacuum contributions to the coefficients in renormalizable 
mean- field models (included in the so-called relativistic Hartree approximation, or RHA [Q]) are 
not natural |41|. In light of our current results in favor of natural coefficients, we view the unnatural 
RHA coefficients as evidence that the vacuum is not treated adequately. 
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to experimental data. At finite density, however, the Fermi momentum enters as a new 
scale. This complicates the situation, because a treatment of the effective lagrangian at 
the level of classical meson fields and valence nucleons (the Dirac-Hartree approximation) 
implies that several different physical effects are implicitly contained in the parameters. 
Not only are short- distance contributions from meson and nucleon vacuum loops absorbed, 
but long-distance, many-body effects from well-known ladder and ring diagrams are in- 
cluded approximately as well |]TB[ . Thus the usefulness of NDA, the naturalness assumption, 
and the expansion in powers of the fields and their derivatives requires further examina- 
tion when these ideas are applied to problems at finite density. It is not obvious that an 
effective lagrangian with natural coefficients that include only short- distance physics will 
remain natural in a one-loop calculation that accounts for higher-order many-body effects 
by adjustments of the coefficients. Here we rely on fits to the properties of finite nuclei to 
determine whether the coefficients are natural and our expansion of the lagrangian is useful. 
(Friar and collaborators have recently investigated naturalness as exhibited by a mean-field 
point-coupling model ^^■) An important topic for future study is to include the many- 
body dynamics explicitly, so as to disentangle it from the short-range and vacuum dynamics 
already implicit in the parameters of the low-energy effective lagrangian. 

The interaction terms of the nucleons and the non-Goldstone bosons can be quite compli- 
cated, and different forms have been used in various models in the literature [53-^. Some 



of these models are actually equivalent if one uses the freedom of field redefinitions [46 



which leave the on-shell S matrix invariant [^7|-|49[|.^ We discuss how to take advantage of 
this freedom to simplify the interaction terms between the nucleons and the non-Goldstone 
bosons. As a result, the lagrangian can be conveniently written in a "standard form", in 
which one mainly has the traditional Yukawa interactions. 

The rest of the paper is organized as follows. In Sec. II, we discuss the construction of an 
effective hadronic lagrangian consistent with chiral symmetry, gauge invariance, and Lorentz 
and parity invariance. Chiral symmetry is realized in the standard nonlinear form of Callan, 
Coleman, Wess, and Zumino (CCWZ) |^8|. Neutral scalar and vector mesons are introduced 



to efficiently describe the mid- and short-range parts of the NN interaction, and the role 
of broken scale invariance in QCD is considered. Georgi's naive dimensional analysis and 
field redefinitions are also discussed, with the new feature of off-shell non-Goldstone bosons. 
In Sec. Ill, we show how the electromagnetic structure of the pion and the nucleon can 
be described in the low-energy regime. Applications to nuclear matter and finite nuclei at 
the Dirac-Hartree level are presented in Sec. IV, and the results are shown in Sec. V. A 
discussion and summary is given in Sec. VI. 



II. CONSTRUCTING THE LAGRANGIAN 

We begin this section by showing how we implement nonlinear chiral symmetry and elec- 
tromagnetic gauge invariance using nucleons, Goldstone pions, and non-Goldstone isovector 
mesons. Next, we discuss the introduction of neutral vector and scalar mesons that are used 



We assume that this is also sufficient for finite-density observables, but we know of no proof. 
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to provide an efficient description of the short- and mid-range nucleon-nucleon interaction. 
To use the resulting lagrangian, which in principle has an infinite number of terms, we must 
have some systematic truncation procedure, and we develop this using the ideas of naive 
dimensional analysis and naturalness, together with some observations on the sizes of meson 
mean fields in nuclei. 



A. Chiral Symmetry 



We illustrate a nonlinear realization of chiral symmetry using a system that consists 
of rho mesons, pions, and nucleons only. We construct a general effective lagrangian that 
respects chiral symmetry, Lorentz invariance, and parity conservation. The role of vector- 
meson dominance (VMD) is also discussed. We restrict ourselves to the chiral limit with 
massless pions and postpone the generalization to a finite pion mass for future study. 

The Goldstone pion fields 7r"(x), with a = 1, 2, and 3 , form an isovector, which can be 
considered as the phase of a chiral rotation of the identity matrix in isospin space: 

^1^ = U{x) =exp{in{x)/f^)lexp{in{x)/f^) . (1) 

Here /jr ~ 93 MeV is the pion-decay constant and the pion field is compactly written as 
7i[x) = tt{x) • |t, with being Pauli matrices. The matrix U is the standard exponential 
representation and the "square root" representation in terms of ^ is particularly convenient 
for including heavy fields in the chiral lagrangian. The isospinor nucleon field is represented 
by a column matrix 




where p{x) and n{x) are the proton and neutron fields respectively. The rho fields p^(a:) 
also form an isovector; we use the notation p^{x) = p^{x) • |t. (See Refs. |]26|j50|j5T| for 



discussions of the formal equivalence between the vector and the tensor formulations for 
spin-1 fields.) 

Following CCWZ we define a nonlinear realization of the chiral group 5'[/(2)l ® 
SU{2)b^ such that, for arbitrary global matrices L G SU{2\ and R G S'f/(2)R, we have the 
mapping 

L®R: (^,p^,iV)— .(e',p;,iV') , (3) 

where 

= Li{x)h\x) = h{x)i{x)R^ , (4) 
p\Xx) = h{x)p^{x)h^{x) , (5) 
N'{x) = h{x)N{x) . (6) 
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As usual, the matrix U transforms as U'{x) = LU{x)K^. The second equahty in Eq. 
defines h{x) imphcitly as a function of L, R, and the local pion fields: h{x) = h{7T{x), L, -R).| 
Under parity, which is defined to be 

V: L^R , L , (7) 

the pion field transforms as a pseudoscalar: ^(t, x) ^''"(t, — x). The existence of the 
parity automorphism (0) allows the simple decomposition U = in Eq. (|l|) [^. It follows 
from Eq. (Q) that h{x) is invariant under parity, that is, h{x) G SU{2)y, with SU{2)y the 
unbroken vector subgroup of SU{2)i^ ® SU{2)^. Eqs. and (|]) ensure that the rho and 
nucleon transform linearly under SU{2)y in accordance with their isospins. Note that the 
matrix h{x) becomes global only for SU{2)y transformations, in which case h = L = R. 

It is straightforward to show that Eqs. @ to @ indeed yield a realization of the chiral 
group [^. We emphasize that there is no local gauge freedom involved, since there are 
no arbitrary local gauge functions in the transformations. The nonlinear realization im- 
plies that chiral rotations involve the absorption and emission of pions, as expected from a 
spontaneously broken symmetry p^|53| . For example, while vector transformations mix the 
proton with the neutron, axial transformations mix the proton with states of nucleons plus 
zero-momentum pions. 

To complete the building blocks of the chiral effective lagrangian, we define an axial 
vector field a^(a;) and a polar vector field v^{x) by 

a,^~{ed,^-^d,e) = al, (8) 

v,^-'^i^%^ + ^d,^^)=vl, (9) 
both of which contain one derivative. Under a chiral transformation, we find 



V 



a 



hv^h^ - ihd^h^ , (10) 
ha^h^ . (11) 



Note that to leading order in derivatives, 

1 i 
= —d^Ti + • • • , v^ = [vr, 9^vr] + ■ • • . (12) 

Jit -^Jtt 

To maintain chiral invariance, instead of using an ordinary derivative acting on a rho or 
nucleon field, we should use a covariant derivative D'^. The chirally covariant derivative of 
the rho field is 



^We can express h in terms of L, R, and U as h{x) = \IU'\x)L\J U{x) = ^ RU^x)L^ L y U{x). 
Given the decomposition L = exp(iQ! • r) exp(i/3 • r), R = exp(— ia • r) exp(i/3 • r), and 
h = exp(i7 • t), with cx, /3, and 7 real, the infinitesimal expansion oi h is -f = (3 — (a x 7r)/2/,r + 
0(«2,/32,7r2) 
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(13) 



We also define the covariant tensors 



Pf,u = DfiPu - DuPf, + ig[Pf,, Pu] ■ 



(14) 
(15) 



The term proportional to [p^, p^] on the right-hand side of Eq. (|TBp vanishes at the mean- field 
level, so the coupling g will be left unspecified. 

To write a general effective lagrangian, we need an organizational scheme for the in- 
teraction terms. We organize the lagrangian in increasing powers of the fields and their 
derivatives. While we do not consistently include loops in this investigation, we assign to 
each interaction term a size of the order of , with k a generic small momentum, and 



d + - + h 



(16) 



where d is the number of derivatives, n the number of nucleon fields, and h the number 
of non-Goldstone boson fields in the interaction term. The first two terms in Eq. (|16D are 
suggested by Weinberg's work |Q. The last term is a generalization that arises because a 
non-Goldstone boson couples to two nucleon fields. Equation (piBf) is also consistent with 



finite-density applications when the density is not too much higher than the nuclear-matter 
equilibrium density, as we will see in Sec. II. D. 

Taking into account chiral symmetry, Lorentz invariance, and parity conservation, we 
may write the lagrangian through quartic order (z/ < 4) as 



£i(x) = N\iY{dfi + + igppf,) + gAl^lbO-^ - M 
-J^Np.^a^'^N-^Nv.^a'^^N 



N 



OTTTT o 



(17) 



where a^i, = z[7^,7j,]/2, g^ ~ 1.26 is the axial coupling constant, M is the nucleon mass, gp 
is the pNN coupling, fp is the so-called tensor coupling, ^fp^vr is the pmr coupling, and P-,,, 
Ktt, h, and I2 are couplings for higher-order vrA^ and tttt interactions. The pion couplings 
can be determined from nN and nn scattering 



T8[] . We have omitted terms with higher 



powers of the rho fields, which have small expectation values in finite nuclei, and four- or 
more-nucleon "contact" terms, which can be represented by appropriate powers of other 
meson fields, as discussed shortly. A third-order {u = 3) term of the form NNii{p^p^) will 
be included by coupling a scalar to two rho fields in Sec. II. E. The antisymmetry of rho 
meson field tensor p^y guarantees three independent degrees of freedom, because there is no 
momentum conjugate to the time component po; cind po satisfies an equation of constraint. 

Chiral symmetry permits isoscalar-scalar and isoscalar- vector terms like (NN)"^ and 
{N^p^Ny. However, analyses of NN scattering with one-boson exchange potentials shows 



8 



that the low-lying mesons are sufficient for describing low-energy many-nucleon systems and 
that they dominate the NN interaction. One may then expect that contact interactions 
of four (or more) fermions can be removed in favor of couplings to massive non-Goldstone 
bosons. The following example shows that this should be a good approximation. 

Consider an additional isovector- vector contact term [A^(t/2)7^A^]^. The contribution 
of this term can be absorbed completely into a renormalization of Qp as follows. If we had 
constructed C by retaining several isovector- vector non-Goldstone bosons, we could eliminate 
them using their field equations to obtain 



where g'^ is the pNN coupling when the heavier mesons are retained, with couplings to 
nucleons given hj gi, ... and masses given by mi , . . . . Evidently, the effects of the heavier 
particles can be absorbed into a redefinition of the pNN coupling: 



9l = 9': + '^9l + --- ■ (19) 



m 



If we approximate the second term in Eq. (113) as 



2 ni p 

the leading error is given by the following term: 
gj mliml - ml) ^ 1 gf 



{-l^){d,mr/2h.N]r , (20) 



9p rnf V 2 



{9,[iV(T/2)7.iV]}^ . (21) 



For a typical resonance in the isovector- vector channel we have mi = 1450 MeV and 



ml{m\ — ml)/m\ \/'o. li gf/g^ is roughly 1/4, the error is about 5% of the already 
small second term. One should be able to test the sensitivity to these errors by allowing the 
non-Goldstone boson masses to vary in the fit to nuclear properties. 
If we were to assume universal VMD in Eq. (0), then 

gpwTT — 9p : ~ ~^ ^ (22) 

and g = gp in Eq. ([T5|) . If we further assume the Kawarabayashi-Suzuki-Riazuddin- 
Fayyazuddin (KSRF) [§^,0 relation 

9p.. - 2^^^, , (23) 

then we may redefine the rho field to be Pp + Vp/gp and the covariant field tensor to be 
P^lu + Vp^/gp, and recover the lagrangian of Weinberg p4[ or Bando et al. [ p^ . 
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B. Gauge Invariance and the Inclusion of the to Meson 



In finite nuclei, electromagnetism also plays an important role. We may include electro- 
magnetic gauge invariance in our formalism in a straightforward manner by noticing how 
the pion, rho, and nucleon fields transform under the local U{1)q symmetry. Here the elec- 
tric charge Q, the third component of the isospin T3, and the hypercharge Y are related by 
Q = T3 + Y/2. 

Under U{1)q, the electromagnetic field transforms in the familiar way 



(24) 



The pion, rho, and nucleon fields transform under the local U{1)q rotation just as in Eqs. 
to (^, with both L and R set equal to 



h{x) 



exp 



Ts + Y 



(25) 



where y = for the pion and the rho, and F = 1 for the nucleon. Indeed, the resulting 
transformations for these hadrons are just as expected from their electric charges. That is. 



[71 , TT ) ^ [IT , e 71 ) 



ipl Pi) - ipl e 



(n, p) — * (n, e p) . 



(26) 
(27) 
(28) 



Electromagnetism breaks chiral symmetry. The lagrangian in Eq. (^^ becomes gauge 
invariant if the ordinary derivatives are replaced with the covariant derivatives 



d^N^[d,+ '-eA^{l + Ts))N 

T-3 



As a result, the axial and vector pion fields become 

v, = v, + \eA,{iq,e]+eq,i]) , 
respectively, and the pion "kinetic" term is 

flir (aX) = f> {aX) - 2ef^A>^tT (v^ts) + 



(29) 

(30) 
(31) 

(32) 
(33) 

(34) 



In addition to terms arising from the preceding minimal substitutions, there can be other 
non-minimal terms in the general effective lagrangian. For example, it is possible to have a 
direct rho-photon coupling term 
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--^F^,tr(r3pn , (35) 

where p'"^ is given by Eq. (^) with the replacements in Eqs. (^) and (0), and F^j^ is the 
usual electromagnetic field tensor. In principle, a similar term proportional to tr {r^W^) 
is also possible, although VMD applied to the pion form factor suggests that the coefficient 
of this term is small Furthermore, the photon can couple to the nucleon as 



~ m^^''^^''"'^ + ^5M^'^'Ar7,(/5, + (3^T,)N + ■ ■ ■ , (36) 

where 

A = ^Ap(l + r3) + ^An(l-r3) , (37) 

with Ap = 1.793 and An = —1.913 the anomalous magnetic moments of the proton and 
the neutron, respectively. The first term in Eq. (^) is needed to reproduce the magnetic 
moments of the nucleons, and the second one contributes to the rms charge radii. Note that 
these terms have z/ = 3 and z/ = 4, respectively. 

We may now write down the resulting manifestly gauge-invariant lagrangian, which con- 
tains the preceding "tilde" fields. Instead of showing this lagrangian in its entirety, however, 
we omit terms that are irrelevant to our applications here. The remaining lagrangian is man- 
ifestly chiral invariant, except for some terms involving the photon field that are proportional 
to the electric charge. We find 

C^{x) = C^{x) - ieA^iV7^(l + r3)iV - -^F^, tr {r-^p^^^) - 2ef^AnT {v,t,) 

-F^J^Xa^'^N + -^5^F'^-iV7,(/5, + /3,r3)iV - -F^.F^^'' . (38) 



AM ^'^ 2M^ ^ - i-y^- ' 4^ 

It is well known that the isoscalar-vector uj meson is needed to describe the short-range 
repulsion of the NN interaction. To have a realistic description of nuclear systems, the 
UJ meson should be included. Note that once the u meson is introduced, an isoscalar- 
vector four-fermion term (A^7^iV)^ can be completely accounted for by adjusting to fit 
experimental data, as discussed above for the rho meson. 

With respect to chiral SU (2) symmetry, the u meson can be treated as a chiral singlet 
represented by a vector field V^{x). The relevant lagrangian can be written as 



^''^''-NV.^a^^N - ^If^^V^" + ■ ■ • , (39) 



AM ^'^ 2g^ 3 

where V^,^ = — 9,^V^. The kinetic terms again refiect the spin-1 nature of the particle. 
Analyses of NN scattering show that the u tensor coupling is not important, since /v is small, 
and this is indeed also the case from our fits (see below). The factor of 1/3 in the coupling 
to the photon arises naturally from SU{3) symmetry and the assumption of ideal mixing of 
the vector (p and lj mesons. The ellipsis represents higher-order terms involving additional 
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powers of the fields and tlieir derivatives, some of which will be considered explicitly in 
Sec. lI.E. 

Before ending this section, we note that it is consistent to integrate out the vector 
0(1020) field because it has a small coupling to the nucleon. The effects of the vector in 
the isoscalar-vector channel can be absorbed into g^, as discussed near the end of Sec. II. A. 
For example, if we integrate out the 0, we obtain a contribution to (3^ of —2^2M'^g^j (Sm^g^) 
and a four-fermion term (N'-y^N)^ with a coefficient proportional to g'^/m'^, where g^f, is its 
coupling constant to the nucleons and is its mass. The contribution of this latter term 
is absorbed into the renormalization of g^. In general, one need only keep the lowest-lying 
meson in a given channel at low energy. 



C. The Isoscalar- Scalar Channel 

Unlike the vector channels, the absence of well-defined, low-lying resonances in the 
isoscalar-scalar channel makes the dynamics more difficult to identify and to model. It is well 
known from the empirical NN scattering amplitude , from one-boson-exchange mod- 



els of the NN interaction ||3^, and from fits to the properties of finite nuclei |^2| that there 



is a strong, mid-range NN attraction arising from the isoscalar-scalar channel, with a range 
corresponding to a mass of roughly 500 MeV. It has also been shown in calculations based on 
either phenomenological scattering amplitudes [^5H37|| or explicit chirally symmetric models 



38| , |39[] that an attraction of the appropriate strength can be generated dynamically by corre- 
lated two-pion exchange between nucleons. No nearby underlying resonance at the relevant 
mass is needed.^ Nevertheless, this dynamics is efficiently, conveniently, and adequately rep- 
resented by the exchange of a light isoscalar-scalar degree of freedom. Therefore, in analogy 
to what we did with the vector mesons, we will introduce an effective scalar meson with 
a Yukawa coupling to the nucleon to incorporate the observed mid-range NN attraction, 
and thereby avoid the explicit computation of multi-loop pion contributions. Based on the 
results of the two-pion-exchange calculations and the success of one-boson-exchange models, 
it is sufficient to give the scalar a well-defined mass, which we will choose here by fitting to 
properties of finite nuclei. Moreover, the light scalar is a chiral singlet, so it does not spoil 
the chiral symmetry of the model. 

In addition to the evidence from the NN interaction, one can ask what other information 
exists on the dynamics in the isoscalar-scalar channel. A natural place to look is the physics 
of broken scale invariance in QCD. This leads to an anomaly in the trace of the QCD 
energy-momentum tensor and also to "low-energy theorems" on connected Green's functions 
constructed from this trace |Q. In pure-glue QCD, these low-energy theorems can be 



satisfied by introducing a scalar glueball field ("gluonium") with a mass of roughly 1.6 GeV. 
When light quarks are added and the chiral symmetry is dynamically broken, the trace 
anomaly can acquire a contribution in addition to that from gluonium. In Ref. |T^, we 
considered this additional contribution to arise from the self-interactions of the light scalar 
mentioned above, and used the low-energy theorems to constrain these interactions. 



^We note, however, that there are some recent indications for a narrow, light scalar |63| 
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The lagrangian, which contained the hght scalar, heavy scalar (glueball), and other 
degrees of freedom, was scale invariant, except for some minimal scale-dependent self- 
interactions among the glueballs and the light scalar. Since the glueballs are heavy, they 
can be integrated out along with other heavy particles and high-momentum modes. If one 
assumes no mixing between the scalar and the glueball, the scalar potential is constrained 
to satisfy the low-energy theorems (at tree level), and the potential is given by [|65| ,p!^ 



Here 5*0 is the vacuum expectation value of the light scalar field S, rris is its mass, and S is 
assumed to have scale dimension d; that is, when x — > X~^x, S{x) X'^S{Xx). 

If we introduce the fluctuation field 0(x) = Sq — S{x) and expand Vs as a power series 
in (j), the coefficients are all determined by the three parameters 5*0, rris, and d. The results 



of Ref. showed, however, that finite-nucleus observables are primarily sensitive only to 
the first three terms in this expansion; thus, even in this extremely simple model of the 
scale breaking, the low-energy theorems do not place significant constraints on the self- 
interactions of the light scalar field. If we allow for the possibility of scalar-glueball mixing 
and employ field redefinitions to simplify the interaction terms between the light scalar and 
the nucleons, the form of the potential Vs can be modified and additional parameters need 
to be introduced. Therefore, to take advantage of a simple Yukawa coupling to the nucleon, 
we adopt here a general potential for the scalar meson, which can be expanded in a Taylor 
series: 

^^ = "^^^2+3! M+4! M^ + ---j- ^^'^ 

Here we have anticipated the discussion of naive dimensional analysis in the next subsection 
by including a factor oil/ f-^ for each power of 0; these factors are then eliminated in favor of 
(7s ~ Af//^. Various counting factors are also included, as discussed in the next subsection. 
If the naturalness assumption is valid, the parameters ^3, ^4, ... should be of order unity. 
This is indeed supported by our fits to finite nuclei, as shown in Sec. V. As with the vector 
mesons discussed earlier, only one low-mass scalar field is needed to a good approximation. 

From a different point of view, chiral symmetry allows for fermion contact terms, such 
as {NNy, [dfj_{NN)]'^, {NNY, and so on. These terms are relevant if there are strong 
interactions between nucleons in the isoscalar-scalar channel. By introducing a scalar field 
with mass mg, Yukawa coupling Qs, and nonlinear couplings K3, K4, ... adjusted to fit 
experimental data, we can account for all of the contact terms. The only relevant questions 
are which implementation allows for the most efficient truncation and how many terms are 
needed in practical calculations. 

We shall not include an isovector-scalar field in our effective lagrangian since the NN 
interaction in that channel is weak p4[. There is no meson with these quantum numbers 
with a mass below 1 GeV, and two (identical) pions in a J = state cannot have T = 1. 
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D. Naive Dimensional Analysis and Field Redefinitions 



The lagrangian constructed thus far is based on symmetry considerations and the desire 
to reproduce reahstic electromagnetic and NN interactions while keeping pion-Ioop calcula- 
tions to a minimum. We have illustrated the lowest-order terms in powers of the fields and 
their derivatives, but in principle, an infinite number of terms is possible, and we must have 
a meaningful way to truncate the effective lagrangian for the theory to have any predictive 
power. 

A naive dimensional analysis (NDA) for assigning a coefficient of the appropriate size to 
any term in an effective lagrangian has been proposed by Georgi and Manohar [ ]5B| , PD[ . The 



basic assumption of "naturalness" is that once the appropriate dimensional scales have been 
extracted using NDA, the remaining overall dimensionless coefficients should all be of order 
unity. For the strong interaction, there are two relevant scales: the pion-decay constant 
fa 93 MeV and a larger scale 0.5 ^ A ^ 1 GeV, which characterizes the mass scale of 
physics beyond Goldstone bosons. The NDA rules for a given term in the lagrangian density 
are 

1. Include a factor of I/Ztt for each strongly interacting field. 

2. Assign an overall factor of /^A^. 

3. Multiply by factors of 1/ A to achieve dimension (mass)^. 

The appropriate mass for A may be the nucleon mass M or a non-Goldstone boson mass; 
the difference is not important for verifying naturalness but can be relevant for N^. counting 
arguments 

As noted by Georgi Rule 1 states that the amplitude for generating any strongly 



interacting particle is proportional to the amplitude /^r for emitting a Goldstone boson. 
Rule 2 can be understood as an overall normalization factor that arises from the standard 
way of writing the mass terms of non-Goldstone bosons (or, in general, the quadratic terms 
in the lagrangian). For example, one may write the mass term of an isoscalar-scalar field 

0(a;) as 

where the scalar mass mg is treated as roughly the same size as A. By applying Rule 1 and 
extracting the overall factor of /^A^, the dimensionless coefficient is of 0(1). Terms with 
gradients (or external fields) will be associated with one power of 1/A for each derivative 
(or field) as a result of integrating out physics above the scale A. (A simple example is 
the expansion at low momentum of a tree-level propagator for a heavy meson of mass mH, 
which leads to terms with powers of d'^ /m\.) It is because of these 1/A suppression factors 
and dimensional analysis that one arrives at Rule 3. 

After the dimensional NDA factors and some appropriate counting factors are extracted 
(such as l/n\ for the n**^ power 0"; see below), the naturalness assumption implies that any 
overall dimensionless coefficients should be of order unity. Without such an assumption. 



14 



an effective lagrangian will not be predictive.[] Until one can derive the effective lagrangian 
from QCD, the naturalness assumption must be checked by fitting to experimental data. 

As an example of NDA and naturalness, we can estimate the sizes of the coupling con- 
stants (?A, Qs, and by applying the rules to the following interaction terms: gj^N'y'^'y^a^N, 
gsN(j)N, and g^N'-j^V^N . Note that the pion field in is already associated with an overall 
factor of I/Ztt [see Eq. ([T^) ] so is simply counted as a derivative. It is straightforward to 
find that 

S-A ~ 1 , g^, g^ r-. ^ < An . (43) 

That is, Eq. (|43|) gives the natural sizes of these couplings, where we use the symbol ~ to 
denote order of magnitude .0 It is reassuring that the preceding values are indeed consistent 
with phenomenology. 

According to NDA, a generic term in the effective lagrangian involving the isoscalar fields 
and the nucleon field can be written as (generalizing to include the pion, rho, and photon is 
straightforward) 

which includes / nucleon bilinears, m scalar fields, n vector fields, and p derivatives. Gamma 
matrices and Lorentz indices have been suppressed. The coupling constant g is dimensionless 
[and of 0(1) if naturalness holds]. 

Numerical counting factors for terms with multiple powers of meson fields are included 
because the NDA rules are actually meant to apply to the tree-level amplitude generated 
by the corresponding vertex. Thus, a term containing 0"*, for example, should have a factor 
1/m! associated with it. These additional factors are certainly relevant when, say, m = 4. 
An explicit illustration of the origin of the counting factors can be found in the pion mass 
term in ChPT: 



\^ir;t. (1/ + l/t - 2) = ^Ifl . (45) 

An important feature of the effective lagrangian is that each derivative on a field is 
associated with a factor of 1/A. Thus organizing the lagrangian in powers of derivatives can 
be a good expansion for describing processes where the dominant energies and momenta are 
not too large. Nevertheless, this organization is intimately connected with the important 



^The assumption of renormalizability would also lead to a finite number of parameters and well- 
defined predictions, but we will not consider this option here [||,|2|. 



^The bound A ^ 47r/,r was originally proposed by Georgi and Manohar |63| and subsequently 



refined to A ^ 47r/^/y']Vy, where Nf is the number of light quark flavors |67|. QCD appears to 
essentially saturate this bound. 



15 



issue of how to remove redundant terms in the effective lagrangian. By employing the 
freedom to redefine the fields, one can organize the lagrangian such that there are no terms 
with (9^ acting on a boson field or 7^9^ acting on a fermion field, except in the kinetic-energy 
terms. These field redefinitions do not change the on-shell S'-matrix elements [^^ -|49[|- 

In finite-density applications, however, the non-Goldstone bosons exchanged by the nu- 
cleons are far off shell, which motivates a different definition of the corresponding field 
variables. Short- distance effects due to heavier particles, high-momentum modes, and dy- 
namical vacuum contributions are assumed to be integrated out, so that their consequences 
appear only indirectly in the coefficients of the lagrangian. Derivatives on the non-Goldstone 
fields are already small compared to their masses in our applications because the momenta 
are limited by the Fermi momentum. Thus, we do not need to eliminate these derivatives. 

Instead, we can use nonlinear field redefinitions to simplify the form of the interaction 
terms between the nucleons and the non-Goldstone bosons so that they take the conven- 
tional Yukawa form. This ensures the consistency of the power counting rule (0). Recall 
that one of the main purposes for introducing the non-Goldstone bosons is to describe the 
intermediate-range NN interactions. Since the nonlinear field redefinitions do not change the 
quadratic kinetic and mass terms, no changes occur to the ranges of the one-boson-exchange 
interactions in free space. All other heavy boson-exchange contributions are short-ranged, 
so we can efficiently parametrize them using meson self-couplings. We can easily show the 
invariance of the one-loop energy functional under field redefinitions, although a general 
proof to all orders seems difficult. It is still convenient to use the field redefinitions to re- 
move all derivatives on the nucleon fields except in the kinetic term. These derivatives can 
be large because valence-nucleon energies are of the order of the nucleon mass M. 

At finite density, the Fermi momentum appears as a new scale and meson fields develop 
expectation values. We must then also decide how to organize the nonderivative terms. 
The key observation is that at normal nuclear-matter density (or in the center of a heavy 
nucleus), the scalar and vector mean-fields are typically 

0.25 M < (?s0o, g.Vo < 0.4 M, (46) 

and the gradients of the mean fields in the nuclear surface, where they are largest, are 

^7s|V0o|, g.\VVo\ ^ {0.2M)^ . (47) 

(Note that the magnitudes of the mean fields are reduced in the surface.) Equations (^Sp 
and (^7]) are our basis for estimating the leading contributions to the energy of any term in 
our lagrangian. 

To summarize, the assumption of naturalness and the observation that the mean fields 
and their derivatives are small compared to M up to moderate densities allows us to organize 
the lagrangian in powers of the fields and their derivatives. We will take advantage of the 
freedom to redefine the fields to simplify the interaction terms between the nucleons and the 
non-Goldstone bosons. For the nucleons, pions, and photons, we use field redefinitions to 
remove high-order derivatives, whereas the scalar and vector meson parts of the lagrangian 
can contain powers of d"^ in addition to the kinetic terms, since these terms are small. By 
truncating the lagrangian at a given order in u, we have a finite number of parameters to 
determine through fitting to nuclear observables. 
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E. The Full Lagrangian 



Based on the preceding discussion, we illustrate here how nonlinear field redefinitions 
allow us to simplify the terms involving interactions between the non-Goldstone bosons and 
the nucleons. Derivatives on the nucleon field other than the kinetic term are removed by 
redefining the field because these derivatives on the nearly on-shell nucleons are of 0{M). 
For the photons and the Goldstone pions, field redefinitions are used to remove higher 
derivatives [^], so that their interactions with the nucleons are unconstrained except by 
gauge invariance and chiral symmetry. Through z/ = 4, the part of the effective lagrangian 
involving nucleons can then be written as 



C^{x) = N{tYV^ + QA^l^a^ -M + g,(l))N 



AM '^^'^ 4M M 

'iu^^^'^^'''^^ - ^^^^^^^ + /5vr3)mF^'^ , (48) 
where the covariant derivative is 

T>^ = df, + iVf, + iQpPf, + ig^Vf, + -eAp{l + rg) . (49) 

One might expect more general isoscalar-scalar couplings to the nucleons. However, all 
possible isoscalar-scalar field combinations of the non-Goldstone boson fields that can couple 
to NN can be redefined in terms of a single isoscalar-scalar field with a Yukawa coupling 
NN(j). For example, the following couplings are redundant: 

NNcj)'^, NN(f)^, NNd'^cj), NNV^V . 

The redefinition of the field induces changes in the coefficients of the mesonic part of the 
lagrangian, but since all terms to a given order in u are already included, there is no loss 
of generality. A similar observation holds for the isoscalar-vector and the isovector-vector 
couplings. 

A term of the form N'-y^^Nd^cj) is irrelevant because the baryon current is conserved. 
Moreover, a coupling of the form 

NNd^V (50) 

can be removed from the hamiltonian by using the constraint equation for the vector field. 
The identity 

9,^^(7,7.5^ + 7.5^7.) (51) 

and the equations of motion allow us to reduce a term of the form NiV^d^N in the la- 
grangian to a sum of Yukawa and tensor terms and the term in ( |50D in the hamiltonian. As 
noted, bilinears in the nucleon fields and derivatives acting on such bilinears can be elim- 
inated in favor of low-mass non-Goldstone bosons. Furthermore, we can use Eq. (^) and 
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partial integrations to reduce mixed-derivative terms such as {Ndfj,N){Nd^N) to the form 
di^{N'j'^'jnN)dcr{N'~f'^'~f^N), which we assume is small compared to terms that are retained 
based on the minor role of tensor mesons in meson-exchange phenomenology. 
In Eq. (^Sp, we have omitted the following z/ = 4 terms: 

Np^.a^^Nfj), NVp^af^'Ncj) . (52) 

The leading effect of the terms in (0) is a slight modification of the tensor couplings fp and 
/v at finite density due to the nonvanishing expectation value of the scalar field. However, 
the effects of the tensor terms are themselves already small, so the further modification is 
negligible. Thus, while we may keep terms through u = 4 in the lagrangian, not all u = 4 
terms contribute equally in a given application. In fact, the pion terms do not contribute to 
the energy of nuclei at tree- level. 

The mesonic part of the lagrangian is also organized in powers of the fields and their 
derivatives. Keeping terms through z/ = 4, we find 

Cu{x) = ^(l + a,^)d,<j)dy + ^ tr (d.Ud^U^) 



e 



2efiAni {v^T,) - —F^, tr {t^p^ + - V^'' 



1 



7 

^2 J,2^ 



+ + + (53, 

where, apart from conventional definitions of some couplings {qs, g^, and g^) and the masses, 
we have defined the parameters so that they are of order unity according to naive dimensional 
analysis. (Again, nuclei will test this hypothesis for us!) Also, since the expectation value 
of the p field is typically an order of magnitude smaller than that of the V field, we have 
retained nonlinear p couplings only through z/ = 3. Note that the ai and a2 terms have 
u = 5. However, from Eqs. ( ^B]) and (^) we can estimate that these two z/ = 5 terms are 
numerically of the same magnitude as the quartic scalar term in the nuclear surface energy, 
so we have retained them. Thus, numerical factors such as l/n\, which are cancelled in 
scattering amplitudes, are relevant in deciding the importance of a term in the energy. 
Our full lagrangian through z/ = 4 is then 

C = C^ + Cm. (54) 

For applications to normal nuclei, keeping these terms should be sufficient if the parameters 
are natural. This may be checked by including the following z/ = 5 terms: 

= ~ 5! '^^ M^^^*^ + 3!^^ ■ 2 ^^^^^ + 4!^^ M ^^^^^^^ " ^^^^ 
As discussed below, the inclusion of £5 improves our fits only marginally. 



18 



III. ELECTROMAGNETIC STRUCTURE OF THE PION AND NUCLEON 



For an effective lagrangian to be useful for calculations of nuclear matter and finite 
nuclei, it should incorporate the low-energy structure of the pion and the nucleon. This 
physics has been extensively studied; there is a long history of vector-meson dominance 
(VMD) descriptions of electromagnetic form factors p8|-pT|. We illustrate here how the 
electromagnetic form factors of the pion and the nucleon at low momenta arise from our 
effective lagrangian through VMD and also through direct couplings of photons to pion and 



nucleon fields [32 



The electromagnetic (EM) current can be obtained from the lagrangian (^), (§5), and 
(^) by taking 5C/5{eA^) after some partial integration: 

+ -(5.Pr + \d,Vn + VMv^r,) . (56) 

Note that the photon can couple to two or more pions either directly or through the exchange 
of a neutral rho meson at non-zero momentum transfers. The coupling of the photon to the 
nucleons is also direct or through the exchange of neutral vector mesons (rho or omega). 

We can determine the tree-level EM form factor of the pion F(— g^) from the current 
(^) and the lagrangian (|5^). For spacelike momentum transfers = —q^, we have 

Thus the mean-square charge radius of the pion is 



dg 

where we have used the experimental values fllGSl 



(0.69 fm)^ , (58) 



2.9, f^ = 2.0, (59) 



47r A-K 



which reproduce the rho widths T p_^t,t, = 150 MeV and FpO^e+e- = 6.8 keV, respectively. 
The observed charge radius is (r 

2)1/2 = (0.66 ± 0.01) fm 
iarly, we can determine th 

nucleon. 



exp 

Similarly, we can determine the EM isoscalar and isovector charge form factors of the 



FnQ') = l-^?^-^7^ + -- ^ (61) 



2 2 M2 2g^ Q2 + ^2 
and the anomalous form factors 
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,2\ \ + fvdv Q 



^(Q^).^-^-^,..., (62) 

+ (63) 

The corresponding mean-square charge radii are 



VM2 3g^ml 



(-%=6f4 + 7^)> (65) 



and 



{r\2 = ^ • (67) 

We will fit and gp to the properties of finite nuclei and then choose /3s, /3v, and /p to 
reproduce the experimental charge radii of the nucleon [pT|: 



(^')sf ~ (^')vf ~0.79fm , (68) 
{r^)l{^ = (0.88 ± 0.05) fm . (69) 

We do not use Eq. (|66D to fix because the remaining radius (r^)g2^ is not sufficiently well 
determined |l3T[| . 



As in Refs. and , we have a contribution from vector dominance and a correction 
to the intrinic structure to second order in a derivative expansion, that is, to order Q^. 
This correction is adequate for our applications to nuclear structure; applications involving 
higher momenta would require us to include additional nonrenormalizable interaction terms 
and loop effects. Thus we do not expect our nucleon form factors to be accurate at large 
momentum transfers. 



IV. FINITE NUCLEI AND NUCLEAR MATTER 



We work at one-baryon-loop order in this paper, which is equivalent to the Dirac-Hartree 
approximation ||6^. This should be sufficient to test the consistency of NDA and the trun- 
cation of our lagrangian. As discussed in Ref. WM, a Hartree calculation can be viewed as 



equivalent to a density-functional approach, in which higher-order many-body corrections 
are treated approximately. Based on Dirac-Brueckner-Hartree-Fock calculations of nuclear 
matter, exchange and correlation corrections to the nucleon self-energies are expected to be 
small The stability of the Hartree results will be tested in a two-loop calculation, which 
will be presented elsewhere. 
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The derivation of the Hartree equations and energy functional for nuclei is similar to 



that described in Ref. [|12|. Because of the large masses of the nucleons and non-Goldstone 
bosons, loop integrals involving these particles include dynamics from distance scales that 
are much shorter than the scale set by the valence-nucleon momenta, which are limited 
by the Fermi momentum kp. As discussed above, these short-range effects are included 
implicitly in the coefficients of the lagrangian. Formally, one can include counterterms to 
remove these loop effects to all orders, which is always possible, since all such terms are 
already contained in the effective lagrangian. 

Our present approach differs from Ref. [0 in that we take a more general point of view 
and employ field redefinitions to simplify the couplings between the nucleon and the non- 
Goldstone bosons. The structure of the nucleon is now described within the theory., so that 
no external electromagnetic form factors are needed, as in earlier calculations [^,^,^]. By 
taking K3, K4, etc. in the scalar potential as free parameters, we also allow for a more general 
description of the dynamics in the isoscalar-scalar channel, as well as for modifications from 
field redefinitions. 



The single-particle Dirac hamiltonian can be written as 
/i(x) = • V + iy(x) + \nR{^) + /?(M - $(x)) + i(l + r3)A(x) 

' ,[3cx • Up\rz^R + /vVTV) + + /^v^s) - -^A/3a • VA , (70) 



where /3 = 70, a = 707, and we have defined scaled mean fields by including the couplings 



W^(x) = g^V,{^) , (71) 

$(x) = ^?s0o(x) , (72) 

i?(x) = gM^) . (73) 

^(x) = eAo(x) . (74) 

The Dirac equation with eigenvalues Ea and eigenfunctions V^a(x) is ||6^,p| 

/iVa(x) = ^a^a(x) , j d^xV-t (x)7/',(x) = 1 . (75) 

We follow the notation and conventions of Ref. [|^ and write the eigenfunctions for 
spherically symmetric nuclei as 

V'a(x) = ^„«;mt(x) = I ^^^"['^l^^"" \Ct, a = {n, K,m}, (76) 



r 



where is a spin spherical harmonic, t = 1/2 for protons, and t = —1/2 for neutrons. 
The radial equations for G and F become 



+ ^)Gair) - [E, - Uiir) + U2{r)]Fa{r) - U,Ga{r) = , (77) 
V,(r) + [E, - f/i(r) - U2{r)]Ga{r) + U,F,{r) = , (78) 



dr r 
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where we have defined the single-particle potentials 

U,{r) = W{r) + taR{r) + (t„ + ^)A{r) + + 2tMV'A{r) , (79) 

U2(r) = M- $(r) , (80) 
Usir) = ^{f^W'ir) + tafpR'ir) + A'{r)[{Xp + A„)/2 + ta{Xp - A„)]} , (81) 

and a prime denotes a radial derivative. 

The mean-field equations for $, W, R, and A are 

-V^Py + mlW = ,^[pB(x) + Av . (p^(x)f)] - + + ^mlW 

^'Jl)^' + ^(V$-VT^ + $V^W^) - ^Pchg(x) , (83) 

-V'R + mji? = i<7?[P3(x) + ■ (pj(x)f)] - r^p^mji? - ^Pchg(x) , (84) 

-V'A = eVchg(x) . (85) 
Here the various densities are defined as 

(86) 
(87) 
(88) 
(89) 
(90) 



Ps(x) 


occ 

= EV^a(x)V'a(x) = 

a 


=f:%^(G?w-f'.V)), 


Pb(x) 


occ 
a 


;x)^a(x) = 


=!:^(G2M+^-.V)), 


Pi(x) 


occ 
a 




*"W = i:^;^;!^2G.(r)F„(r), 


P3(x) 


occ 
a 


;x)r3^a(x; 


) = !:^(2M(GSW + f.^W), 

a " ' 


pj(x) 


occ 
a 


;^x)ir3/3a • 


fV^„(x) = ^^Z^(2t„)2G„(r)F„(r) 



where the summation superscript "occ" means that the sum runs only over occupied (va- 
lence) states in the Fermi sea. Evidently, these densities depend only on the radial coordi- 
nate. The charge density is given by 
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Pchg(x) = Pd(x) + Pm(x) , (91) 

where the "direct" nucleon charge density is 

Pd(x) = pp(x) + • (pl(x)f) + ^[/JsV^pB + /^vV^s] , (92) 

and the vector meson contribution is 

p„,(x) = -^V^R + -^V^W . (93) 

Here the "point" proton density and nucleon tensor density are determined (as always) from 
the solutions to the Dirac equation in the mean fields: 

OCC 

Pp(x) = -E^IW(l + ^3)^a(x) = -(PB + PS) , (94) 

OCC 

pj(x) = Y.ijl{^)t\pcx.fM^) , (95) 

a 

with the anomalous magnetic moment A from Eq. (0). As seen in the next section, the 
vector meson contribution and the final term in Eq. ([9^ ) generate the nucleon form factor 
that reduces the fiuctuations of the point-nucleon charge densities. Note that the mean-field 
equations can be used to eliminate the second derivatives in the charge density in favor of 
first derivatives, which are more accurately calculated numerically. 
The energy functional is given by 

OCC „ 

E=Y.E^- / d'xU^ , (96) 

where 

f/m = -\^Ps + \W{PB + pI) + ^i?(p3 + pj) + \Apa 

%2 ]vf Vi2 ^ 24M 80M3/ 4 M ^2 



4Mi^^ + ^^M + yM^^^^ 

— fCo + ^-W^ + ^-{V^f - ^-{VWf (97) 



and we have applied the mean-field equations to arrive at this form. 

The center-of-mass (cm.) correction to the nuclear binding energy can be estimated 
nonrelativistically using 

- 2MS ' ^^^^ 
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where B = Z + N is the baryon number of the nucleus, and 



(^cm) = -E(«lv» + EK«lv|/?)r , (99) 

a a,f3 



with 



{a\d\f3) = J d^xilji{^)dM^) ■ (100) 

For simphcity, we use an empirical estimate given by Reinhard |[70| : 

Ecu = (MeV) , (101) 

which determines (-Pcm) ^^^^ from Eq. (p^) . The binding energy e is then 

e = BM-{E- Ecu) ■ (102) 

Note that we cannot replace BM by Zrrip + Nrrin, the sum of the masses of the protons and 
neutrons, unless we use different masses for the proton and neutron in the Dirac equation. 

There are also cm. corrections to the nuclear charge density and charge radius. We 
estimate these in analogy to the well-known nonrelativistic results for harmonic-oscillator 
wavef unctions, which are exact. The charge density in momentum space, i.e., the charge 
form factor, becomes 

i^chg(q) = ^Pchg(q)(l + + ■■■), (103) 

where Pchg(q) is the Fourier transform of the charge density pchg(x). We expect this cor- 
rection to be accurate for momenta |q| ^ 2fm~^. With the cm. correction, the predicted 
charge density in coordinate space is, to a good approximation, 

{pchg(x)}obs ~ Pchg(x) - A VVchg(x) . (104) 
The mean-square charge radius is given by 

3 



(r^)chg = (r^) - , (105) 



where 



(r') = |/d3xxVchg(x). (106) 

The energy density for uniform, symmetric nuclear matter in the one-loop approximation 
can be obtained from the preceding results by observing that the single-particle energy 
eigenvalue becomes 



E{k) = W + \/k^ + M*^ , (107) 
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where M* = M — $, and $ and W are now constant mean fields [see Eqs. ([72|) and ([fl|)]. 
The energy density £ becomes 



£[M*,p3] = -^$^(- + -- + -- + 



gl V2 ' 3! M 4! M2 5! 

- 4!^(^° + ^4)^' + + (2^ /''^'^ V'^^T^ , (108) 

where fcp is the Fermi momentum and pe = ^k'^/Sir'^ is the baryon density. An extensive 
analysis of energy densities of this form is given in Ref . |T^ . One can also compute the bulk 
symmetry-energy coefficient [jl[ 



«4 = TT^j-ofcg + Z / . ^ (109) 



0^ 1 F 



where the effective rho mass m* is defined here by 

m!' = mj(l + r]p^/M) . (110) 



V. RESULTS 



To test the physical ideas discussed in the preceding sections, we fit our one-loop results 
to the properties of finite nuclei. We choose a set of spherical nuclei and solve the Dirac 
and mean-field equations for each nucleus by iteration, until self-consistency is reached ^2 



We then calculate a set of observables {X^lj} for each nucleus and tune the parameters to 
minimize the generalized lEIl defined by 



X 



EE 

i X 



cxp 



(111) 



where i runs over the set of nuclei, X runs over the set of observables, the subscript "exp" 
indicates the experimental value of the observable, and are the relative weights. The 
nuclei we choose are ^"^O, ''"Ca, ^^Ca, ^^Sr, and 2°^Pb. 

Similar fitting procedures have been standard for determining parameter sets for rela- 
tivistic mean-field models [^,^,^,0. In the previous studies, however, the principal goal 
was to determine minimal parameter sets that provided a good description of bulk nuclear 
properties, so as to maximize predictability. The set of observables used to make the fit was 
also kept at a minimum. In these fits, the parameters included only the scalar mass mg, the 
nucleon-meson couplings Qs, g^, and Qp, and two scalar self-interaction parameters, K3 and 
K4. (The exception is the point coupling model of Ref. ||71|, which has a corresponding set 



of nucleon-only point couplings.) The nucleon and vector meson masses were generally fixed 
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at experimental values. The resulting parameter sets provide good fits to nuclear binding 
energies, charge distributions, and single-particle structure, for a wide variety of nuclei. 

Our effective lagrangian, truncated at u = 4, contains many additional parameters. Thus 
we expect to find good fits to nuclear properties. However, our goals are quite different from 
the previous studies. We are not concerned with maximizing predictive power; instead, we 
want to use experimental data to test whether the parameters of the effective lagrangian are, 
in fact, natural. In so doing we want to gain insight into an appropriate level of truncation. 
Therefore we constrain the parameters by using every observable we believe should be well 
reproduced. 

The function in Eq. ( |111| ) serves as a figure of merit for our fit. We do not expect 
our model to reproduce the experimental observables at the level of the actual experimental 
uncertainties in each case; this means that a direct statistical interpretation of will not 
be meaningful. Rather we choose the weights to be the relative accuracy we expect in a 
good fit. A "perfect fit" corresponds to a contribution of unity to from each observable. 
In choosing weights we follow closely the prescription in Ref. so that we will have a 



basis of comparision. Different choices of weights will emphasize different aspects of the 
physics and lead to different "optimal" parametrizations. However, we have tried a range of 
reasonable weights and the qualitative conclusions given below are independent. 

The choice of observables starts with the basic physics of the semi-empirical mass formula, 
which should be reproduced by any reasonable model of nuclei. Here this is achieved directly 
by fitting to the binding energies of individual nuclei. We adopt a relative weight of 0.15% 
for each binding energy. This implies that binding energies per nucleon should be reproduced 
to two decimal places, with an uncertainty of roughly unity in the second place. But it is 
also important that the systematics of the energies are followed; otherwise, the minimization 
procedure may find reasonable fits to energies but with poor systematics (particularly when 
fewer parameters are used, as with = 2 or = 3 truncations). Therefore we also use the 
deviation of experimental and calculated binding energies 6ei = (ei)exp — (ei)th to determine 
6ai, 6a2, and Sa^ according to 

6ei = 6ai Ai - da^ A]''' - Sa^ (N^ - Z^f/Ai . (112) 

Here Ni and Zi are the number of neutrons and protons in the i^^ nucleus and Ai = Ni + Zi. 
Thus an exact fit to the energies would also have 5ai = 6a2 = Sa^ = 0. The deviations 6a2 
and 5a4 are included as separate terms in Eq. ( |111| ) in the form [^o-i/M/^aJ^, with Wsai = 0.08. 
(This weight was selected based on the results for low-order truncations; it has no intrinsic 
significance.) 

We also expect the effective theory to accurately reproduce low-momentum observables, 
which leads us to focus on the charge form factor. We choose the rms charge radii of the 
nuclei and the so-called diffraction- minimum-sharp (d.m.s.) radii as observables (another 



choice would have been the surface thickness |j73[). The d.m.s. radius of a nucleus is defined 
to be |7| 

i?dms = , (113) 

where Qo^^ is the three-momentum transfer at the first zero of the nuclear charge form 
factor F{Q) = Fchg(q) with Q = |q|. The d.m.s. radius is expected to be a clean low- 
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momentum quantity for a nucleus and should be well described. The weights chosen for the 
radii correspond to an accuracy of roughly two decimal places. 

Finally, we expect to reproduce various features of the single-particle structure in nuclei. 
We focus in particular on spin-orbit splittings near the Fermi surface, which are reasonably 
well known for both proton and neutron levels. These splittings are closely related to the 
size of the mean fields at the one-loop level. In a previous study, a good reproduction of 



these splittings was found even when they were not included as observables [|73|. However, 
our expanded parameter space means that this correlation may not follow without explicitly 
including them. In addition, we use the proton I/19/2 energy level in ^osp]^ ^^le overall 

scale of the energy levels; the positioning of this level helps to determine parameters that 
influence the symmetry energy. Finally, the proton splitting Ep{2d3/2) — Ep{lhii/2) is chosen 
as a characteristic splitting between states with large and small /, which constrains the shape 
of the effective single-particle potential |7^. The additional weight on observables for lead 



also leads to a more sensitive determination of the isovector parameters. 

To summarize, the observables (a total of 29) and their relative weights are taken as 
follows: 

• The binding energies per nucleon e/B, with a relative weight of 0.15% 

• The rms charge radii {r'^)chg^ with a relative weight of 0.2% 

• The d.m.s. radii -Rdms, with a relative weight of 0.15% 

• The spin-orbit splittings AEso of the least-bound protons and neutrons, with a relative 
weight of 5% for ^^O, 15% for ^ospb, 25% for ^^Ca and ^^Ca, and 50% for ^^Sr 

• The proton energy Ep{lhQ/2) and the proton level splitting Ep{2d3/2) — -Ep(l/iii/2) in 
208pb, with relative weights of 5% and 25%, respectively 

• The surface-energy and symmetry-energy deviation coefficients 5a2 and Sa^ [see 



Eq. ([ml) ], each with a weight of 0.08. 



For further discussion of the relative weig hts for e/B, (r^)^^ , and AEgo, see Ref. [fn] 
Additional observables that might be added in future investigations include neutron radii, 
quadrupole deformations, monopole resonance energies, and an enlarged set of test nuclei. 

We take the nucleon, uj, and p masses to be given by their experimental values: M = 
939 MeV, rriy = 782 MeV, and rrip = 770 MeV.0 The anomalous magnetic moments of the 
nucleon are fixed at Ap = 1.793 and An = —1.913. We also fix = 5.01 as given in Eq. (p^. 
The free-space charge radii of the nucleon, as given in Eqs. (^) and (0), are used to fix Ps, 
/3v, and fp by solving Eqs. (|64D to (|67D. The remaining thirteen parameters Qs, g^, gp, ?7i, 



^Includina; free parameter starting from the parameter sets in Table | leads to only a tiny 

improvement in (less than unity) and a small change in my (about 10 MeV). Therefore we fix 
the oj and p masses, since the latter also has little effect. 
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TABLE I. Parameter sets from fits to finite nuclei. 





W\ 


C\ 


Ql 


Q2 


Gl 


G2 


rUs/M 


0.60305 


0.53874 


0.53735 


0.54268 


0.53963 


0.55410 




0.93797 


0.77756 


0.81024 


0.78661 


0.78532 


0.83522 


5v/47r 


1.13652 


0.98486 


1.02125 


0.97202 


0.96512 


1.01560 


ffp/47r 


0.77787 


0.65053 


0.70261 


0.68096 


0.69844 


0.75467 


m 




0.29577 






0.07060 


0.64992 


m 










-0.96161 


0.10975 


«3 




1.6698 


1.6582 


1.7424 


2.2067 


3.2467 








-6.6045 


-8.4836 


-10.090 


0.63152 


Co 








-1.7750 


3.5249 


2.6416 


'Hp 










-0.2722 


0.3901 


Oil 










1.8549 


1.7234 


0.2 










1.7880 


-1.5798 


/v/4 










0.1079 


0.1734 


/p/4 


0.9332 


1.1159 


1.0332 


1.0660 


1.0393 


0.9619 


/3s 


-0.38482 


-0.01915 


-0.10689 


0.01181 


0.02844 


-0.09328 




-0.54618 


-0.07120 


-0.26545 


-0.18470 


-0.24992 


-0.45964 



7^2, Tip^ K3, K4, Co, '"^s; /v; c^i; and ^2 for the V = Al parametrization are then obtained by the 
optimization procedures described above. 

This optimization program faces serious difficulties. With such a large parameter set 
and because the observables themselves are highly correlated, there is a definite problem of 
underdetermination. The minimization routines used here are reliable for determining local 
minima but are not guaranteed to find a global minimum, so we cannot be sure to have 
sampled all relevant values of the parameters. There are further complications in obtaining 
a more global fit because the required stability of the Hartree iterations precludes large steps 
through parameter space. Thus the minimization routines must be efficient in navigating 
narrow valleys in parameter space. To make our results as robust as possible, we used two 
different codes and made many runs with different starting parameters. 

Our studies of the full z/ = 4 parametrization yielded two distinct parameter sets (Gl 
and G2) with essentially the same values.[^ These sets and several truncated parameter 
sets described below are presented in Table |. The observables of the five nuclei for the 
Gl and G2 parameter sets are summarized in the tables and figures. Bulk binding-energy 
systematics are given in Table The rms charge radii and d.m.s. radii are summarized in 
Table we have consistently used the experimental Fourier-Bessel charge densities and 
charge radii of Ref. for all nuclei. The spin-orbit splittings of the least-bound nucleons 
are given in Table For comparison, we also include results from the point-coupling 



model of Ref. |71]. Note that our rms radii contains cm. corrections, while those defined in 



Ref. [71 1 do not. The charge densities (and point-proton densities) and charge form factors 



Small changes in the weights can make the lower for one or the other, but not significantly. 
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TABLE II. Binding-energy systematics for sets Gl, G2, and the point-coupling model of 
Ref. ||7l| (set PC). Binding energies per nucleon are given in MeV. 



Set 40ca 

~G1 7^98 8!55" 

G2 7.97 8.55 

PC 7.97 8.58 

expt. 7.98 8.55 



48 Ca 88Sr 208pb 

8.67 8.72 7.88 

8.68 8.72 7.87 

8.69 8.75 7.86 
8.67 8.73 7.87 



TABLE III. Rms charge radii (in fm) and d.m.s. radii (in fm) for sets Gl, G2, and the 
point-coupling model of Ref. (set PC). 



Set 


160 




40 


Ca 


48 


Ca 


88Sr 




208 


Pb 






Rdms 




Rdms 




Rdms 








-Kdms 


Gl 
G2 
PC 

expt. 


2.72 
2.73 
2.73 
2.74 


2.77 
2.77 

2.76 


3.46 
3.46 
3.45 
3.45 


3.84 
3.84 

3.85 


3.45 
3.45 
3.48 
3.45 


3.95 
3.95 

3.96 


4.19 
4.19 
4.21 
4.20 


4.99 
5.00 

4.99 


5.50 
5.50 
5.50 
5.50 


6.80 
6.80 

6.78 



for 1^0, 40Ca, ^SQa, and ^ospb are shown in Figs. I to| 
and 208p]^ g^j^g presented in Figs. |^ to 11. 



The single-particle spectra for Ca 



The fits to nuclear charge radii, d.m.s radii, binding energies, and spin-orbit splittings 
are quite good; they are almost all at the relative accuracy prescribed by the corresponding 
weights. A good reproduction of the spin-orbit force in finite nuclei necessarily leads to large 
scalar and vector mean fields in the interiors of nuclei or in nuclear matter near equilibrium 
density. Our calculations show that our parameters produce scalar and vector mean fields 
consistent with the estimates given in Eqs. (^) and (^Tj). The single-particle spectrum in 
lead is also quite good, although the level density near the Fermi surface is somewhat too 
low, as usual in the Dirac-Hartree approximation. 

The low-momentum behavior of the nuclear charge form factors is well reproduced, while 



TABLE IV. Spin-orbit splittings for the least-bound neutrons and protons for sets Gl, G2, and 
the point-coupling model of Ref. (set PC). 







neutron AE^so (MeV) 






proton A^so (MeV) 






Gl 


G2 


PC 


expt. 


Gl 


G2 


PC 


expt. 


160 


6.0 


6.0 


6.4 


6.2 


6.0 


6.0 


6.4 


6.3 


40 Ca 


6.6 


6.5 


6.8 


6.3 


6.6 


6.5 


6.8 


7.2 


48Ca 


5.8 


5.6 


5.9 


3.6 


6.3 


6.1 


6.1 


4.3 


88Sr 


2.1 


2.0 


1.9 


1.5'^ 


6.2 


6.0 


6.1 


3.5" 


208pb 


0.8 


0.9 


0.9 


0.9 


1.8 


1.8 


2.0 


1.3 



"Calculated from Bohr and Mottelson ||7^ as in Ref. [|7T 
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FIG. 1. Charge density of ^^O. The sohd hne is taken from experiment ||7^. Charge densities 
are shown for the Gl and G2 parameter sets from Table I. Also shown is the point-proton density 
[pp of Eq. for set Gl. 
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FIG. 3. Charge density of ^''Ca. The sohd hne is taken from experiment |77]. Charge densities 
are shown for the Gl and G2 parameter sets from Table I. Also shown is the point-proton density 
for set Gl. 
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FIG. 5. Charge density of ^'^Ca. The sohd hne is taken from experiment |77]. Charge densities 
are shown for the Gl and G2 parameter sets from Table I. Also shown is the point-proton density 
for set G2. 
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FIG. 7. Charge density of '^°Pb. The sohd hne is taken from experiment [77|. Charge densities 
are shown for the Gl and G2 parameter sets from Table I. Also shown is the point-proton density 
for set G2. 
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FIG. 9. Predicted single-particle spectra for Ca using the Gl and G2 parameter sets from 



Table I. The leftmost values are from experiment |7£]. 



the high-momentum parts show some departure from the experimental data. This is ex- 
pected, since our approximations are not as accurate at higher momentum transfers and 
we have only constrained the low-momentum form factor. In coordinate space, the vec- 
tor meson and gradient contributions to the charge densities are important in reducing the 
point-nucleon charge fluctuations. This is evident in the figures, where the point-proton 
density is also shown. We emphasize that the charge density is obtained directly from the 
Hartree solutions; there is no convolution with external single-nucleon form factors. 

Experience with a broad class of relativistic mean-field models shows that models that 
successfully reproduce bulk properties of finite nuclei share characteristic properties in infi- 
nite nuclear matter |Tl],|13|. These properties are the equilibrium binding energy and density, 



TABLE V. Nuclear matter equilibrium properties for sets from Table | and for the 
point-coupling model of Ref. ||7l| (set PC). Values are given for the binding energy per nucleon 
(in MeV), the Fermi momentum kp (in fm~^), the compression modulus K (in MeV), the bulk 
symmetry energy coefficient 04 (in MeV), M* /M, and QvVo (in MeV) at equilibrium. 



Set E/B-M kp K 04 M*/M q^Vq 

Wi -16.46 1.279 569 409 0.532 363" 

CI -16.19 1.293 304 32.0 0.657 255 

Ql -16.10 1.299 242 36.4 0.597 306 

Q2 -16.13 1.303 279 35.2 0.614 292 

Gl -16.14 1.314 215 38.5 0.634 274 

G2 -16.07 1.315 215 36.4 0.664 248 

PC -16.13 1.299 264 37.0 0.575 322 
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FIG. 10. Predicted neutron single-particle spectra for Pb near the Fermi surface using the 
Gl and G2 parameter sets from Table I. The leftmost values are from experiment |7S], and PC 



stands for the point-coupling model of Ref. 1 71 ] 



the compression modulus K, the value of M*/M at equilibrium, and the symmetry energy 
04. One further condition is that the light scalar mass satisfies 500 ^ rris ^ 550 MeV. This 
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FIG. 11. Predicted proton single-particle spectra for ^'^'^Pb near the Fermi surface using the 



Gl and G2 parameter sets from Table I. The leftmost values are from experiment |7S], and PC 
stands for the point-coupling model of Ref. |71|. Note that the I/19/2 level is the lowest level in the 
first unoccupied shell. 
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FIG. 12. Contributions to the energy per particle in nuclear matter for parameter sets Gl and 
G2 from the n*'^-order terms of the form ^'■W"^, where n = l + m. The boxes are terms with / = 0, 
the circles are terms with m = 0, and absolute values are shown. Results from set Gl are open 



and those from G2 are filled. The crosses are estimates based on Eq. (|44|). The arrow indicates 
the total binding energy eo = 16.1 MeV. 



condition ensures reasonably smooth charge densities, good surface-energy systematics, and 
an appropriate shape for the single-particle potentials [Tl|]. This is confirmed from the prop- 



erties of nuclear matter predicted from our fits to finite nuclei, as shown in Table ^ for the 
parameter sets in Table |. 

It is instructive to examine the contributions to the energy of nuclear matter from various 



terms in the energy density. This is shown in Fig. |T2|, where the contributions to the energy 
per nucleon at equilibrium from the quadratic, cubic, and quartic powers of the fields $ and 
W in Eq. ( |108| ) are indicated separately. We set the scale for the natural size of contributions 
expected from an n*'^-order term using the power counting estimate [Eq. (Q)] for an n^^- 
order contact term and the approximation (NN) ^ {N'^N) = p^, with pa being the baryon 
density and the bracket indicating an expectation value. These estimates are marked with 
crosses, and the error bars reflect a range of A from 500 MeV to 1 GeV. We observe that 
the estimates are consistent with the field energies found from the fits. Furthermore, these 
estimates also set the scales for the mean meson fields, which are consistent with Eq. (^6|). 
This is an indication that the naturalness assumption is valid (see below). 
Several observations can be made based on the figure: 

• The contributions decrease steadily with increasing powers of the fields, so that a 
truncation is justified. Note that this decrease will become more gradual as the density 
is increased above equilibrium density. 

• The nuclear matter binding energy (indicated by the arrow) is an order of magnitude 
smaller than the individual u = 2 contributions. This implies that the u = 2 scalar 
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and vector contributions cancel almost completely and also that z/ = 4 contributions 
are still quite significant (at least for set Gl). 



The z/ = 5 contributions [see Eq. (|53|)] are representative results obtained when the 
= 4 fit is extended. (They correspond to the lowest values of x^)- They are not 
well determined; almost identical values of can be found for very different z/ = 5 
parameters. Nevertheless, even when unnatural z/ = 5 parameters are introduced, their 
net contribution to the energy is typically a few tenths of an MeV. Thus a truncation 
at z/ = 4 is consistent and adequate. 

The systematics of the two parameter sets at z/ = 4 are quite different. Set Gl parallels 
the systematics of the standard parameter sets from Refs. Il73|j70| , [7^ , while the z/ = 4 



terms for G2 are quite small. This is possible because of the contributions from the ai 
and a2 terms, which involve gradients of the meson fields and thus do not contribute in 
nuclear matter. The contributions of these gradients [note the signs of the and the 
signs in Eq. (^)] to the spin-orbit sphttings in finite nuclei allow set G2 to accurately 
reproduce the empirical spectra, even though M*/M in nuclear matter is larger than 
in conventional models. 

If we return now to the parameters in Table |, we see that while the Yukawa couplings and 
the scalar mass are close for the two sets, the meson self-couplings are very different. This 
is consistent with our earlier Hartree calculations []62| , p!3| , which showed that the bulk and 



single-particle observables of interest are determined to a large extent by a relatively small 
number of nuclear properties. As noted above, by reproducing five important properties of 
nuclear matter and by choosing an appropriate scalar mass, one essentially fixes the predicted 
properties of doubly magic nuclei, and thus the thirteen parameters in the present model are 
underdetermined. The parameters are therefore highly correlated and are sensitive to the 
relative weights chosen for the observables, including the spin-orbit splittings. To determine 
the parameters better, one must include additional observables that are sensitive to different 
aspects of the dynamics. 

To explore these issues further and to put our z/ = 4 fit in perspective, it is useful to 
find parameter sets for more constrained truncations. If we optimize with the "standard" 
mean- field model, which includes the Yukawa couplings, the scalar mass, and and K4, 
we find ~ HO; ^ind the parameter set (which we label Ql) is qualitatively close to the 
"best-fit" sets of the past and also to the corresponding parameters of set Gl (see Table |). 
Adding the z/ = 4 vector self-interaction, with parameter (q, improves the to about 90 
(set Q2). The x"^ values of sets Gl and G2 are slightly less than 50, so it is clear that a 
better fit is obtained with more parameters, but not dramatically so. 

If we truncate our model at z/ = 2 (set Wl), keeping only the Yukawa couplings and 
the scalar mass and setting all other parameters to zero, the best is over 1700. The 
"optimal" set in this case depends strongly on the relative weights, because only a subset 
of the observables can be reproduced accurately. The standard choice [^ has been to 
reproduce the charge radii accurately at the expense of significant underbinding. The present 
set of weights leads to much better binding energies (despite a large compressibility) but the 
systematics of the charge radii (oxygen in particular) is badly reproduced. The spin-orbit 
splittings are too large in all cases, corresponding to a low value for M*/M (see Table |Vp. 
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At z/ = 3 (set CI), we add two parameters, K3 and 771. (We hold r]p at zero for simplicity.) 
If the number of parameters was the sole key to a good fit, we would expect to be 
comparable to that from the standard models. However, the best fit has ~ 400. In 
practice the rji contribution plays the same role as K4 if one looks at energy contributions in 
nuclear matter, but it is insufficient to accurately reproduce nuclear sizes. It is likely that 
the density dependence of the two = 3 terms is too similar to provide enough flexibility 
for a good fit. 

Thus we find that we need at least z/ = 4 terms to find a good fit to nuclear prop- 
erties. Examples are sets Ql and Q2. Once we allow all u = 4 terms, however, the fit 
is under deter mined, as is evident by the differences between our sets Gl and G2. One 
of the important new features are the ai and ^2 parameters, which accompany derivative 
corrections. Without these included in the fit, set G2 cannot achieve a better than 100. 
If we include z/ = 5 parameters from Eq. ( ^5]) in the fit, only small improvements in 
(< 2 units) are found. The net contribution to the energy from these terms is very 
small (roughly a few tenths of an MeV) in these fits, even if the coefficients are artificially 
made large so that individual contributions are significant. Furthermore, the changes in the 
lower-order coefficients from the u = 4 fits are small. Thus we find z/ = 4 contributions 
are necessary and sufficient for this set of observables. Clearly the u = 5 results, when 
included, do not drive the physics, and their effects can be absorbed into small changes in 
the lower-order couplings. The insensitivity to the u = 5 terms validates our truncation of 
the lagrangian. 

The parameters in Table | have been displayed in such a way that they should all be of 
order unity according to NDA and the naturalness assumption. This is seen to be the case.|^ 
It is important to emphasize this nontrivial result: Without the guidance from NDA, these 
coefficients could be any size at all! We conclude that NDA and the naturalness assumption 
are indeed compatible with and implied by the observed properties of finite nuclei, even 
though we are absorbing many-body effects in the coefficients. 

However, the detailed mechanism by which the nuclear observables constrain the param- 
eters to be natural is not transparent. Second-order [u = 2) coefficients vary only slightly, 
regardless of the number of parameters used. The other parameters, while natural, are not 
as well determined. We searched for unnatural parameter sets by starting from unnatural 
values, but the optimization always produced natural results. The z/ = 5 parameters can be 
made large, but they are tuned by the optimization so that they cancel almost completely 
in the Hartree equations. Despite the narrow range of densities present in ordinary nuclei, 
the different density dependence of various terms in the energy is apparently sufficient to 
enforce naturalness. 



The most unnatural coefficient is K4 from set Gl. Its contribution to the energy is still compatible 
with natural estimates despite its size because of the 1/4! counting factor. 
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VI. DISCUSSION AND SUMMARY 



Our goals in this paper are to construct an effective field theory based on hadrons that 
is appropriate for calculations of finite-density nuclear systems and to propose methods for 
making these calculations tractable. The lagrangian is consistent with the underlying sym- 
metries of QCD, and heavy, non-Goldstone bosons are introduced to avoid the calculation of 
(some) virtual pion loops. To organize the lagrangian, we rely on Naive Dimensional Analy- 
sis (NDA) , which allows us to extract the dimensional scales of any term, on the assumption 
of naturalness, which says that the remaining dimensionless coefficient for each term should 
be of order unity, and on the observation that the ratios of mean meson fields (and their 
gradients) to the nucleon mass (which generically represents the "heavy" mass scale) are 
good expansion parameters. Since the meson fields are roughly proportional to the nuclear 
density, and since the spatial variations in nuclei are determined by the momentum distri- 
butions of the valence nucleon wave functions, this organizational scheme is essentially an 
expansion in kp/M, where kp is a Fermi wavenumber corresponding to normal nuclear den- 
sities. If NDA and the naturalness assumption are valid, one can expand the lagrangian in 
powers of the fields and their derivatives, truncate at some finite order, and thus have some 
predictive power for the properties of nuclei. Our fits to bulk and single-particle nuclear 
properties show that this is indeed the case. 

As in all effective field theories, the coefficients in our nonrenormalizable lagrangian im- 
plicitly include contributions from energy scales heavier than our generic mass A. Thus 
contributions from baryon vacuum loops and non-Goldstone boson loops (including, in prin- 
ciple, particles more massive than those considered here) are parametrized in this way. 
Moreover, since we presently compute the nuclear observables at the one-baryon-loop or 
Hartree level, fitting the parameters to nuclei means that they also implicitly include the 
effects of many-body corrections. These consist of nucleon exchange and correlation con- 
tributions, which may involve virtual states with large excitation energies, modifications of 
the vacuum dynamics at finite density, like the nuclear Lamb shift, and virtual pion loops, 
since the mean pion field vanishes in spherical nuclei. These many-body effects involve both 
short- and long-range dynamics, so it is notable that our fitted parameters are still natural. 

We can understand this result with the following arguments. First, it is unlikely that 
natural parameters arise from a sensitive cancellation between unnatural parameters in the 
lagrangian and compensating large many-body effects. It is much more likely that the 
parameters in the lagrangian are natural, and the modifications from many-body effects are 
small. This is consistent with explicit calculations of relativistic exchange and correlation 
effects, which show that these contributions do not significantly modify the nucleon self- 
energies nor introduce large state dependence, at least for states in the Fermi sea.Q This 
"Hartree dominance" of the mean fields [jl3| and of the parameters needed to produce them 



^^A well-known, but somewhat exceptional example concerns the symmetry energy in nuclei. 
Explicit calculations show that one-pion exchange graphs make a significant contribution to the 
symmetry energy |8C], but these effects can be simulated at the Hartree level by simply choosing 
a rho-nucleon coupling that is larger than what one would expect from free-space considerations. 
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is a crucial element in the expansion and truncation scheme described above. An important 
topic for future work is the explicit evaluation of correlation effects in this model, both to 
disentangle these contributions from the fit parameters and to verify the extent to which 
Hartree dominance is valid. 

We can also understand why earlier relativistic Hartree approaches to nuclear structure 
were STicccssfTil. Many of these calculations used lagrangians that were truncated arbitrarily 
at quartic interactions in the meson fields, but we have seen that quintic and higher-order 
terms are unimportant (and undetermined) for the nuclear observables of interest. Moreover, 
even though most of these investigations arbitrarily set some of the cubic and quartic meson 
parameters to zero, we have also seen that the full complement of parameters through fourth 
order (u < 4) are underdetermined by the data. Therefore, keeping only a subset of the 
parameters does not preclude the possibility of a realistic fit to nuclei. 

In summary, we have constructed an effective lagrangian that maintains the symmetries 
of QCD, such as Lorentz invariance, parity conservation, chiral symmetry, and electromag- 
netic gauge invariance. The bulk and single-particle properties of nuclei can be described 
realistically by this effective hadronic theory at the one-baryon-loop level. Naive dimensional 
analysis and the naturalness assumption are found to be compatible with the observed nu- 
clear properties. The fitted parameters are natural, and the results are not driven by the 
last terms kept. The description of the spin-orbit splittings of the least bound nucleons and 
the shell structure are generally good, although the errors are relatively large compared with 
the other observables. The long-range electromagnetic structure of the nucleon is included 
explicitly using vector-meson dominance and derivative couplings to the photon; it is not 
necessary to introduce ad hoc electromagnetic form factors. The mesonic and gradient con- 
tributions are important for reducing fluctuations in the point-nucleon charge density, and 
the resulting nuclear charge form factors agree well with experiment. The extension of these 
calculations beyond one-baryon-loop order is currently in progress. 
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